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, Abstract 
^ ■ 

^ I We show how the hypothesis that galaxy discs conform to self-similar dy- 

,— , ■ namics leads to the identification of an annular region of the optical disc 

, which is such that, corresponding to the TuUy-Fisher scaling law defined 

^ \ on the exterior annular boundary, there is a similar scaling law defined on 

^ ' the interior annular boundary. This result is confirmed at the level of sta- 

Q^. tistical certainty over several large ORC samples. Furthermore, the same 

Q \ analysis provides insight into the uncertainties associated with the "best" 

^ ' way of defining Vmax, the rotation velocity used for the Tully-Fisher scaling 

^ ■ law and Rmax, the galaxy radius at which Vmax is measured. Finally, as a di- 

I— I \ rect consequence, we are led to a refined Tully-Fisher law which is largely 

^ I insensitive to the means by which a galaxy's rotation velocity is defined. 

> ■ 

^ : 

J2 : 1 Introduction 

p ; 

00 ■ 1.1 General comments 
o ; 

^ ■ One of the many problems which hinders our understanding of the underly- 

^ ■ ing primary processes driving spiral galaxies is the fact that such galaxies are 

• rH . frequently observed to be "afflicted" by one of the following: 

><: 

. 1. ongoing interactions with external objects; 

2. manifest signatures of such interactions in the near past; 

3. internal inhomogeities generating local perturbations; 

4. presence of bars; 
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5. unusually active central regions, 

and so on. The net effect of these various phenomena is to considerably com- 
plexify the task of identifying the irreducible physics and phenomenology which 
define the essential nature of the spiral galaxy. 

1.2 Structure suggesting self-similar dynamics in annular discs 

Regardless of the foregoing problems, there still exist various (statistically) strong 
signatures which appear to be fundamental; these are: 

1. the existence of the Tully-Fisher scaling relationship (M ^ aQ + bo log Vmax)r 

2. the long-standing recognition (Danver 1942, Kennicutt 1981) that the spi- 
rality of the spiral arms in disc galaxies can be usefully classified in terms 
of the logarithmic spiral, R = Roexp{b6), for radial displacement R and 
angular displacement 6. 

The second of these signatures, in particular, is consistent with the hypothesis 
that the dynamics in spiral discs conform to laws of self-similarity. That is, in the 
annular sub-regions where spirality is manifest, 

yrot yrot , , 
= -. — = ... = constant, 

yrad yrad 

from which, immediately, 

Vrot = Af{R), Vrad = Bf{R) (1) 

for some undetermined function f{R) and constants (A, B). There are some 
obvious comments to be made: 

1. Since spirality is a property restricted to an annular region of galaxy discs, 
and since the (logarithmic) spirality property is implied by dH), then this 
latter hypothesized velocity field must be considered confined to some 
annular region. The precise quantitative definition of this annular region 
is given in §3. It 

2. There is a large amount of galaxy rotation data available which makes a 
statistical study of y^f perfectly feasible so that, in principle, it is possible 
to hypothesize and test specific proposals for f{R); 
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3. Until recently, it was thought that there were essentially no systematic 
radial flows in galaxy discs - presumably because of the mass-flow im- 
plications - but it is becoming increasingly recognized that radial flows 
are a common feature of galaxy discs so that, in principle at least, we can 
conceive the possibility Vj-^d <^ Vwt for the typical galaxy disc. However, 
the very small amount of radial flow data presently available renders any 
large scale statistical study of radial flows impossible, for some years at 
least. 

1.3 An hypothesis for f{R) and the manner of its testing 

The topic of how best to express rotation curves in some generic functional form 
has received much attention over the years; but all of these efforts have been fo- 
cussed on the whole rotation curve. By contrast, here we are concerned only with 
the rotation curve in a limited annular region (to be precisely defined in §3.1|) 
of the optical domain where self-similar behaviour appears to be manifested - 
that is, we are not concerned with rotation curve behaviour in the interior re- 
gions, and nor are we concerned with the behaviour where they tend to become 
(approximately) flat. The most simple possibility that has any chance of accom- 
modating the variation between discs for the two-dimensional velocity field in 
the optical annulus only is f{R) = i^*, so that our hypothesis becomes 

Vrot ~ AR«, V,,d ~ BR'', R^,n <R< Rmax 

where Rmin ^rid Rmax defined the boundaries of the annular region to be de- 
fined in §3.11 However, as we have already indicated, it is only feasible to con- 
sider whether or not Vyot ~ AJ^* at the moment. We assess the validity of the 
power-law hypothesis in qualitative and quantitative analyses and briefly de- 
scribe these in the immediately following. 

1.3.1 Qualitative assessment 

The adoption of the power-law hypothesis Vrot = -^R"^ for the annular disc leads 
directly to the following results: 

1. the classical Tully-Fisher scaling relationship, 

Mtf ^ao + bo log Vmax, 

defined on the exterior boundary of the annular disc, emerges automati- 
cally from the correlation which is found to exist between the power-law 
parameters {A, a.) and global luminosity properties; 
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2. the power-law exponent, a, which can be estimated independently of dis- 
tance scales being determined, allows the definition of an enhanced Tully- 
Fisher relation of the type 

Mtf = a + (& + ac) log Ymax / 

where Vmax can be defined in many ways without appreciably affecting 
the quality of the resulting M^f determinations; 

3. corresponding to the TuUy-Fisher scaling relationship on the exterior bound- 
ary of the annular disc, there is a directly analogous scaling relationship 

Mff ^ fll + &1 log Vynin 

defined on the interior boundary of the annular disc. This result is entirely 
new, and lead to the identification of the dynamically coherent annular disc 
as an objective reality. 

The fact that these quantitative results arise directly from the hypothesis of 
Vrot = AR* for the annular disc provides considerable circumstantial support 
for the hypothesis. 



1.3.2 Quantitative assessment 

Generally speaking, when considering power-law fits to data, it is not possible 
to say anything more than the power-law does /does not provide a good fit. In 
the present case, for various reasons, we can go much further: we are able to 
say that, for all practical purposes, the data behaves as if Vrot — AR* is the fun- 
damental law and not merely just a good fit. 

This is done in the following way: we note that any whole-disc rotation curve 
can be expressed as 



vb V^oy " \RoJ 

for arbitrary values {Rq > 0, Vq > 0, a.) and for some suitably chosen function, 
g{R/Ro). We are able to show that for each ORG in the large samples consid- 
ered, values of (Rq, Vq, a) can be chosen (uniquely for each disc) such that 



log 



N(x,cr) where x ^ 0, a <<1, 



over the annular discs of the galaxies concerned. In other words, g{R/Ro) itself 
behaves similarly to a log-normal random variable with a mean close to zero. 



4 



For example, for the case of the Mathewson, Ford and Bucchom (1992) sample 
of 900 ORCs, we fmd N{x,(r) ^ N(0.002, 0.07). 

The implication of this latter result is that, over large samples, the rotational 
dynamics on the annular disc (to be given an objective operational definition) 
are described to high precision in a least-square sense by V ^ [^]modei = 
where (A, a) are parameters determined by the ORC in the annular disc. Al- 
ternatively, we can say that deviations from simple power law behaviour in the 
annular disc are normally distributed with a mean close to zero. 

2 The rotation curve samples 

The optical rotation curve samples considered in this analysis are those of: 

1. Mathewson, Ford & Bucchom (1992) which consists of 900 objects here- 
after referred to as MFB; 

2. Mathewson and Ford (1996) which consists of 1200 objects, hereafter re- 
ferred to as MF; 

3. Courteaux (1997) which consists of 305 objects, hereafter referred to as SC; 

4. A composite sample collected from Dale, Giovanelli and Haynes (1997, 
1998, 1999) and Dale & Uson (2000) consisting of 497 objects, hereafter re- 
ferred to as DGHU. This sample was provided by permission of Giovanelli 
and Haynes. 

2.1 The folding of rotation curves 

The analysis to be described is predicated upon the ability to fold rotation curves 
with high precision and, in this, Persic, Salucci & Stel (Persic, M., Salucci, P., & 
Stel, F 1996, MNRAS, 281, 27) laid out the basic approach to be adopted. 

Traditionally, the 'quick and dirty' approach to folding involves simply iden- 
tifying the photometric centre of a disc (an unambiguous process) and then 
subtracting the observed redshift at that point from the whole rotation curve. 
This approach is perfectly acceptable if one is primarily interested in the bulk- 
flow of galaxies. But, if ones primary interest is in the internal dynamics of the 
galaxies themselves (as was Persic et al's), then it is necessary to recognize that 
the centre of rotation of a galaxy (its kinematic centre) does not necessarily co- 
incide with its photometric centre. Persic et al found that the best results were 



5 



obtained if the centre of folding (the kinematic centre) and the systematic reces- 
sion velocity were treated as two free parameters, to be determined such that 
the 'folded arms were maximally matched'. 

The difficulty, of course, is in deciding what is meant by 'maximally matched' 
- especially given that rotation curves are, more often than not, unevenly sam- 
pled in some way; for example, the measurements may not extend equally on 
either side of the disc, or measurements may be taken at unequal intervals on ei- 
ther side of the disc, etc. Persic et al adopted a labour-intensive, non-automated 
and qualitative 'eye-ball' approach to the problem, whilst Catinella, Haynes and 
Giovanelli (AJ, 130, 10371048, 2005) have constructed a method based upon us- 
ing the two parameters to optimize the fit of Persic et al's Universal Rotation 
Curve function to their rotation curves. 

The automatic algorithm developed by this author (A&A Supp, 140, 247-260, 
1999) is also based upon Persic et al's two-parameter method and, briefly, oper- 
ates as follows: 

The two parameters are varied such that, in a Fourier decomposition 
of the whole rotation curve (which is purely asymmetric in the ideal 
case), a normalized measure of the residual symmetric components 
is minimized. 

Due to various problems, such as asymmetric sampling or how many Fourier 
modes to use for any given ORG etc, the actual implementation of this simple 
algorithm is complicated, but the code is available from this author. 

Finally, one very important point, identified by Persic et al, is that, given a whole 
set of velocity measurements over a rotation curve it is necessary to also have 
quantitative estimates of the absolute accuracy of each individual velocity mea- 
surement available - and then to use this information as a means of filtering out 
only the best individual measurements for the folding process. Broadly speak- 
ing, any individual velocity measurement is retained only if its estimated abso- 
lute error < 5%. For the MFB, MF, SC and DGHU samples, this requirement 
led to losses of 35%, 25%, 46% and 46% respectively of all individual velocity 
measurements. The net effect of these data losses meant that many ORCs were 
then left with insufficient data points on them to permit a reliable folding. The 
overall attrition rate of ORCs lost to the overall analysis via this process were 
3%, 4%, 7% and 8% respectively 
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3 The annular disc: A dynamically coherent whole 



The most significant result of this section is the demonstration that the annu- 
lar disc is a dynamically coherent and objectively defined component of the to- 
tal disc. Apart from giving the algorithmic definition of the annular disc, the 
demonstration consists in showing how the interior boundary of the annular 
disc is a dynamical transition region between the inner and outer parts of the 
disc on which a scaling law, directly analogous to the classical TuUy-Fisher scal- 
ing law, is defined. 

3.1 Extraction of the annular disc from the whole disc 

The process to be described is based on the hypothesis that V f« [V]^^^^^ = AR* 
describes rotation velocities in ORCs in certain annular regions which are to be 
extracted from the whole disc. Before describing the extraction algorithm, there 
are two important notes to be made: 

1. By the 'whole disc' in this context, we mean the complete set of veloc- 
ity data for the disc with the exception of any filtered-out poor-quality 
individual measurements (cf last section). The annular disc is extracted 
from the whole disc using a statistical algorithm based upon the technique 
of linear regression: following conventional definitions, an observation is 
reckoned to be unusual if the predictor is unusual, or if the response is un- 
usual. For a p-parameter model, a predictor is commonly defined to be 
unusual if its leverage > 3p/N, when there are N observations. In the 
present case, we have a two-parameter model so that p = 2. Similarly, the 
response is commonly defined to be unusual if its standardized residual 
>2. 

2. The prescription V ^ AR* can only be completely specified for any given 
ORC when the distance scale has been set. However, the exponent a. is 
independent of the distance scale; consequently, the identification of any 
annular region over which the approximate power-law may be valid is 
also independent of the distance scale. For this reason, the extraction al- 
gorithm described below assumes that radial displacements are specified 
in radians, R(radians) say; the corresponding calculated A is denoted as A*. 

The computation of (a. In A*), for any given folded and inclination-corrected 
rotation curve, can now be described using the following algorithm.: 

1. Assume the model V ~ \y\ynodei ~ ^*^lradians) rotation velocity in 

some annular region of the disc, to be determined. But initially, use the 
whole disc ORC; 
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2. Form an estimate of the parameter-pair (a. In A*) by regressing the InV 
data on the lnR(^„^,„„s) data for the folded ORC; 

3. Determine if the innermost observation only is an unusual observation in 
the sense defined above; 

4. If the innermost observation is unusual, then exclude it from the compu- 
tation and repeat the process from (2) above on the reduced data-set; 

5. If the innermost observation is not unusual, then no further computation 
is required: 

(a) the remaining set of points at this stage are considered to define the 
extracted annular region. The minimum value of R for any given 

disc, Rffiin 

say, defines the radius of the interior boundary of the ex- 
tracted annulus and V^in denotes the rotation velocity on this interior 
boundary; 

(b) the current value of a defines the exponent of the power-law over the 
extracted annular region; 

(c) the value A* has no particular significance since a distance scale has 
yet to be set. 

This algorithm has the result that, on average, (a. In A*) is computed on the 
exterior 88% of the data points in each ORC of the MFB sample, the exterior 87% 
of data points in each ORC of the MF sample, the exterior 91% of data points in 
each ORC of the DGHU sample and the exterior 91% of the data points in each 
ORC of the SC sample. 

3.2 An old result quantified and refreshed 

It has been known for a long time that, in a qualitative sense, the steeper the rise 
of a rotation curve on its interior portion, then the flatter the rotation curve is in 
its exterior portion. 

Given our hypothesis that V ~ AR'^ on the annular disc, then it is clear that 
the parameter A has the numerical value of the approximate rotation velocity at 
Ikpc (given that R is in units of kpc) - as estimated by using the power law to 
extrapolate (or interpolate) the data to Ikpc. It follows that the value of A can 
also be considered as a proxy measure for the steepness of the initial rise of the 
rotation curve in which case, given the qualitative statement above, we should 
expect A and a to be in an inverse relation to each other. 
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MFB (1992) data. 
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Figure 1: 



For the sake of demonstrating the reality of this relationship, we use MFB's 
Tully-Fisher distances to set the radial distance scale so that A can be computed 
for each galaxy in sample. Figure [T] then shows the scatter diagram plotting 
(a. In A), computed as in §3.11 for each foldable galaxy in the Mathewson et al 
(1992) sample. The powerful nature of their inverse relation is clearly apparent. 

To summarize, the given analysis quantifies old qualitative knowledge by pro- 
viding an explicit and well-defined partition of the optical disc into distinct dy- 
namical regions: the interior disc containing the steep initial rise, and the exte- 
rior annular disc over which - as we shall demonstrate - the rotational dynamics 
are described to extremely high fidelity by the V ~ AVS^ . 

3.3 Scaling laws and the annular disc 

The algorithm described in §3. II will always produce a result in the form of an 
extracted annular disc. So, the significant questions are: 

1. Does the process described have any basis in physics so that the annular 
disc can be identified as a physically coherent distinct sub-component of 
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the whole disc ? 

2. Or is the data reduction process essentially arbitrary, the product of which 
has no physical significance whatsoever? 

In the following, we show that a scaling law similar to the TuUy-Fisher law 
which applies at the exterior boundary of the annular disc also applies on its in- 
terior boundary, R = Rmin determined in §3.1[ and at similar levels of statistical 
significance. Thus, this interior boundary is a boundary of objective physical 
significance which scales according to luminosity and which defines, in effect, 
the transition within the disc between one form of behaviour and another. 



3.3.1 The interior Tully-Fisher law 

In the following, we establish the existence of a TuUy-Fisher-type law on the 
interior boundary of the annular disc. The results of the individual regressions 
are given below: 



MFB: Mtf ~ (-14.60 ± 0.15) + (-1.40 ± 0.03) In V^/„, 



n = 849, ^^,„d = -42, =67%; 



MF: Mtf ~ (-14.79 ± 0.13) + (-1.36 ± 0.03) In V^/„, 

^adj 



n = 1070, tgrad = -45, = 65%; 



(3) 

SC: Mtf ~ (-14.89 ± 0.26) + (-1.36 ± 0.05) In V,«,„ , 



n = 275, t,rad = -25, 69%; 



DGHU: Mtf ^ (-11.81 ±0.24) + (-1.95 ±0.05) In y„„„, 
n = 487, tgrad = -3S, Rlaj = 75%. 

The quoted i-statistic is each case refers to the estimated gradient value and it is 
quite obvious that, in every case, the dependence of the Tully-Fisher estimated 
absolute luminosity, Mtf> upon the rotational velocity on the interior boundary 
of the annular disc is established at the level of statistical certainty. 



It is interesting to note the remarkable consistency in the regression models 
arising from the MFB, MF and SC samples, with that arising from the DGHU 
sample being significantly different, the reasons for the discrepant nature of the 
DGHU results are not clear. 
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3.4 Conclusions for the annular disc 



It has been demonstrated, at the level of statistical certainty, that a scaling law 
analogous to that of TuUy-Fisher, exists on the interior boundary on the ex- 
tracted annular disc. This establishes that: 

1. the extraction process is not arbitrary, but gives a result which has an ob- 
jective physical significance; 

2. the annular disc is a physically coherent distinct component of the com- 
plete disc; 

3. the boundary between the inner-disc and the annular disc is a physical 
transition boundary between one form of dynamical behaviour and an- 



4. the power-law hypothesis, V = AR*, for the annular disc provides a good 
statistical description for the rotational dynamics in annular discs over 
large samples. This conclusion will be quantified in a later section. 

4 An enhanced Tully-Fisher law: the parameter oc 

The classical Tully-Fisher scaling law is given by 



where Vrot is an estimated maximal (or characteristic) rotation velocity. The ac- 
tual precise definition of Vj-ot is highly problematical - especially when working 
in the optical - and very much effort has been expended in trying to construct al- 
gorithmic definitions for it. In this section, we consider an enhanced Tully-Fisher 
scaling law 



where the parameter a is the exponent in the power-law fit Vrot = ^R'^ to any 
given rotation curve, and we show that the inclusion of the term a log Vyot in dl]): 

• makes the performance of the Tully-Fisher scaling law insensitive to the 
definition of Vrot', 

• improves the performance of the Tully-Fisher scaling law when the chosen 
definition of Vrot is sub-optimal. 



other; 



Mtf = a + blog, Vrot 




(4) 
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4.1 Vrot definitions used in the present analysis 

The various definitions of Vrot used in this analysis to make the point are listed 
as follows: 

• The histogram method, for which Vyot = Vhist = 0-5 (V^ — ^ioo-n) where 
Vj^ is the N-th percentile velocity. Versions of this are used by MFB, MF 
and DGHU; 

• The model linewidth method, for which Vyot = Vmodei where Vmodei is the in- 
terpolated velocity at some formally defined position on some phenomeno- 
logically defined function fitted to the rotation curve. SC investigated a 
variety of these to arrive at his Vmax (which is the peak velocity on his 
model curve) and his V2.2 (which is his model velocity at K = 2.15 x (disc 
scale length)). Catinella et al (2005) use their polyex function to arrive at a 
similar definition for Vmodei- 

• Various authors (eg Persic & Salucci 1995) have suggested using Vrot = 
Vopt where Vopt is the (usually interpolated) rotation velocity at the optical 
radius, defined as Rgpt = ^83- In the following analysis, Ropt{radians) is 
determined directly from the photometry and we estimate Vopt by fitting 
Vrot = ^*^radians ORG data prior to the distance scale being set (a. is 
independent of the distance scale) and then defining Vopt = ^*^1pt(radians) ■ 

4.2 The analysis 

For each of the four samples, MFB, MF, DGHU and SG, available to us we: 

• calibrate the classical Tully-Fisher scaling law using Hubble distances to- 
gether with: 

- the rotational velocities provided by the authors concerned. For MFB, 
MF and DGHU, these are estimated using versions of the histogram 
method so that Vrot = Vy^ist- For SG, these are estimated using various 
model linewidths; we consider just two of these, his V2.2 and Vmax', 

- the rotation velocities interpolated to the optical radii, RgS/ ^opt say; 

• calibrate the enhanced Tully-Fisher scaling law, dD, for the same samples 
and definitions of Vrot', 

• compare the fits of the classical Tully-Fisher scaling laws to the enhanced 
scaling laws over each of the four samples and for each of the definitions 
of Vrot considered. 
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Efficiency of classical TF vs enhanced TF 



MFB 

Vhist Vopt 



DGHU 

Vhist Vopt 



MF 

Vhist Vopt 



sc 

V2.2 Vmax Vopt 



# = Classical Tully-Fisher 
= Enhanced Tully-Fisher 



Figure 2: 



The detailed analyses are given in the appendices but the results are summa- 
rized in figure |2] : The samples are given in the major columns and the major 
columns are subdivided according to which definition of Vyot is used. The rel- 
ative scaling-law performance is estimated using a normalized measure of the 
fit of the scaling law to the data, as follows: 

• for each sample, the classical scaling law is calibrated using V^/g; in the 
case of MFB, MF and DGHU and V2.2 in the case of SC; 

• in each case, an estimate of fit to the data is given by the adjusted 
parameter, labelled Rj^pg, Rj^iP' -^dghu ^^'^ ^sc respectively; 

• each of R%ipB' ^mf' ^dghu ^^'^ ^sc normalized to unity; 

• for every other calibration of either the classical scaling law or the en- 
hanced scaling law, the corresponding estimate of fit, R^, is normalized 
using one of R%ipp, R%ip> ^dghu ^sc appropriate. 

In figure |2l the filled circles respresent the normalized -R^ values arising from 
the classical scaling law whilst the open diamonds represent these values aris- 
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ing from the enhanced scaling law. The figure makes it very clear that the effec- 
tiveness of the enhanced scaling-law is very much less sensitive to the definition 
of Vrot than that of the classical scaling law. In particular: 

• in every case, Vrot = Vopt performs poorly in the classical scaling law but 
performs comparably with all other choices in the enhanced scaling law; 

• for the DGHU and SC samples, the choices V^fgf and ^2.2 respectively ap- 
pear to be optimal in the sense that the use of the enhanced scaling law 
makes no measurable difference over the classical scaling law; 

• for the MFB and MF samples, the enhanced scaling law gives measurable, 
but small, improvements over the classical scaling law when Vj^igt is used. 
The detailed analysis, given in the appendices, shows that these improve- 
ments are strongly significant in a statistical sense; 

• for the SC sample, the enhanced scaling law gives much larger measurable 
improvement over the classical scaling law when V2.2 is used. 

4.3 Interim conclusions 

We have demonstrated that the inclusion of the predictor a log Vrot has the po- 
tential to greatly improve the predictive power of the TuUy-Fisher scaling law 
but that, in practice, improvements are strongly contingent of the precise def- 
inition of Vrot employed. It appears that there is an optimal definition of Vrot 
which is similar to those employed by DGHU and SC but, whichever definition 
is used, suboptimal choices are brought to near-optimality by the inclusion of 
a log Vrot 3S a predictor. 

5 The classical TuUy-Fisher scaling law 

So far, we have shown that new insights into scaling laws for galaxy discs arise 
when we assume that V ^ AR* over some partition of galaxy discs, R^in < 
R < Rmax say. In the following, we: 

1. show how the TuUy-Fisher scaling law emerges from the power-law hy- 
pothesis for the annular disc; 

2. clarify why a In Vrot can be expected to be a significant predictor according 
to the precise definition chosen for Vrot- 
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5.1 The power-law hypthesis and the TuUy-Fisher scaling law 



Given the hypothesis V = AR'^ for the rotation velocity in the annular disc, 
where (A, a.) vary between discs, then for an arbitrarily chosen point {Rq, Vq) 
on any given rotation curve, we have 

InA = InVo - alnRo- (5) 

But we know that if, for example, we choose Vq = Vrot on any given ORG then 
correlations of the type Vq = Vyot = /(M, a) exist. 

In practice, an exploration of the MFB, MF, DGHU and SG samples (see ^for 
the quantitative numerical details) shows that models of the general type 

In A ^ [lnA]^^^g; = uq + a^M + a {h + biM + bshiS) , (6) 

where S is the surface brightness, are comprehensive (that is, all significant pre- 
dictors are included) and account for, typically, 95% of the variation in In A data. 

A quick comparison of (jSj with Q shows that the decomposition 

lnVo = flo + «iM, InRo = - (^1 + + fcalnS) (7) 

is possible so that, immediately, we have the classical TuUy-Fisher scaling law, 
together with a corresponding relation giving the value oiR = Rq at which Vq is 
measured. However, the decomposition ^ is not unique and in this observation 
we find an understanding of the appearance of a. In Vrot as a predictor. 

5.2 The predictor a In Vrot in the Tully-Fisher scaling law 

The considerations of §5.11 above lead to an understanding of why a In Vrot is 
also a predictor in Tully-Fisher calibrations, and of why its significance varies 
according to the definition of Vrot employed. 

In effect, (O defines both a particular characteristic velocity on an ORG and the 
characteristic radius at which this velocity is to be measured, both being deter- 
mined in terms of luminosity properties of the galaxy concerned. However, 
suppose that we wish to defined the characteristic velocity at 

InRo = lnRo + AlnRo 

rather than at In Rq- Then (O shows that, since M and S are fixed for the galaxy, 
AlnRo 7^ can only arise from some perturbation of (&i,&2/&3) - for example, 
A In Rq = — A&2M. In this case, (jS) can only be satisfied if (jZll becomes 

In Vq = «o + («i - ^^Afc2) M, In Rq = - {h + + Ab2M + &3 InS) . 
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Thus, the simple act of redefining the characteristic radius according to Rq Rq 
immediately causes a; In Vq to become a predictor for M. 

In other words, whether or not a. In Vq is a significant predictor in any given case 
depends entirely upon the definition adopted for Vq which, in effect, amounts 
to systematically selecting a particular definition for Rq. 

6 The detailed analysis of the four samples 

In the following, we show how the considerations of the previous section work 
in practice in the context of the samples of MFB, MF, SC and DGHU: in particu- 
lar, we explore the model 

In A ^ [In A]^„^g; = uq + uiMtf + a (&i + biMTF + &3 In Stf) (8) 

for those samples where, in each case, Mjp is estimated using Tully-Fisher rela- 
tionships calibrated by the authors concerned. 

6.1 The sample of MFB 

The MFB sample contains 864 ORCs that we were able to fold and, after remov- 
ing 28 of the more extreme outliers, the model © is defined by the table: 



Coeffs 


ttQ a\ h\ 1)2 ^3 




Estimate 
Std Error 
i-statistic 


-0.946 -0.286 5.590 0.475 0.586 
0.149 0.007 0.286 0.015 0.020 
-6 -41 20 31 29 


95% 



so that, from (O we immediately get 

In Vo = -0.946 - 0.286MrF, In J^o = -5.590 - QA7bMTF - 0.586 In Sjf. (9) 
The first of these two relations gives, directly (to two decimal places) 

Mtf = -8.051ogio Vq - 3.31 
which is very close to MFB's calibration of the Tully-Fisher relation, 

Mtf = -7.961ogio Vrot - 3.30. 
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6.2 The sample of MF 



The MF sample contains 1083 ORCs that we were able to fold and, after remov- 
ing 26 of the more extreme outliers, the model dH]) is defined by the table: 



Coeffs 


uq a\ h\ b2 ^3 




Estimate 
Std Error 
f-statistic 


-1.199 -0.298 5.976 0.469 0.477 
0.141 0.007 0.288 0.015 0.019 
-9 -45 21 32 25 


94% 



so that, from (O we immediately get 

In^o = -1.199 - 0.298Mtf, InJ^o = -5.976 - 0.469Mtf - 0.4771nSTf. (10) 

The first of these two relations gives, directly (to two decimal places) 

Mtf = -7.45 log^o Vq - 4.62 

which is to be compared with MF's (MFB's) calibration of the TuUy-Fisher rela- 
tion, 

Mtf = -7.961ogio Vrot - 3.30. 
6.3 The sample of DGHU 

The DGHU sample contains 497 ORCs that we were able to fold and, after re- 
moving 24 of the more extreme outliers, the model (H)) is defined by the table: 



Coeffs 


flO ^\ bi &2 ^3 




Estimate 
Std Error 
f-statistic 


-1.412 -0.307 5.540 0.434 0.415 
0.180 0.008 0.420 0.021 0.024 
-8 -37 13 20 17 


96% 



so that, from ^ we immediately get 

In Vo = -1.412 - 0.307Mtf, InRo = -5.540 - 0.434Mtf - 0.415 In Stf- (H) 
The first of these two relations gives, directly (to two decimal places) 

Mtf = -7.50 log^Q Vq - 4.60 
which is to be compared with DGHU's calibration of the Tully-Fisher relation, 

Mtf = -7.68 logio^^of- 4.23. 
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6.4 The sample of SC using Vmax 

The SC sample contains 282 ORCs that we were able to fold and, after removing 
24 of the more extreme outliers, the model ^ is defined by the table: 



Coeffs 


aQ a\ hi b2 




Estimate 
Std Error 
f-statistic 


-3.225 -0.394 9.282 0.604 0.288 
0.219 0.010 0.744 0.0235 0.026 
-15 -39 12 17 11 


98% 



so that, from (O we immediately get 

In Vo = -3.225 - 0.394MrF, Ini^o = -9.282 - 0.604Mtf - 0.288 In Srf. (12) 
The first of these two relations gives, directly (to two decimal places) 

Mtf = -5.841ogio Vo - 8.19 
which is to be compared with SC's calibration of the TuUy-Fisher relation using 

his Vffiax/ 

Mtf = -e.l9\og^QVrot -7.50. 

7 A quantitative test of the power-law hypothesis 

So far, the assumption of the power-law hypothesis, Vrot = AR"^, has led to con- 
siderable new insight into the scaling properties of galaxy discs so that ipso facto 
the hypothesis has great utility in practice; however, these results amount to 
qualitative evidence in favour of the hypothesis. 

In the following we give a quantitative data analysis which shows that, at the 
level of statistical certainty, rotation velocities within the annular disc (defined 
in @ of spirals behave as if the power law hypothesis 

~ 1 / ^min ^ R ^ Rfnax / 
0/ 

where {a, Vq, Rq) are parameters unique to each galaxy, is the physical law gov- 
erning rotation velocity within spiral discs. 
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Sample size = 453 



- Best fit normal: Log(X) ~ N(0.00 1,0.033) 
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Best fit normal: Log(X) ~ N(0.002,0.026) 
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7.1 The analysis 

Any rotation curve can be expressed as 



V 



(13) 



for arbitrary (Rq > 0, Vq > 0, a.) and a suitably chosen function, g (R/Rq). From 
this there follows immediately 



V 



R_y 



R 



MS 



RoJ 



log 



V 



a log +logX. 



In the following, we show that, on the basis of the four large samples, X 
a statistical sense. The analysis proceeds as follows: 



(14) 



1 in 



choose the scaling functions relating {Vq,Rq) to luminosity properties, as 
one of (HJ, (|TOl), dH]) or ((121) depending on the ORC sample being consid- 
ered; 

with the chosen scaling functions, perform a linear regression of log (V/Vq) 
on log {R/Rq) as described in §3.11 



we now have a linear model for (|T4|) for the chosen ORC where the gradi- 
ent gives an estimate of a and the zero point represents the log X term. 

the normalized histograms of the log X distributions, together with the 
best fitting pdf, is given in figure 2. For each of the four samples, the best 
pdf fits are given in the table: 



Sample 


N{x,a) fitted to logX 


MFB 


N(-0.001, 0.043) 


MF 


N(+0.001, 0.039) 


SC 


N(+0.002, 0.026) 


DGHU 


N(+0.001, 0.033) 



• It follows from the latter table and ((14)) that X = 1 in the statistical sense, 
and at the level of near statistical certainty. In other words, there is abso- 
lutely no evidence to suggest that the power-law hypothesis for rotation 
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velocities in the annular disc is not supported and that, for all practical 
purposes, we have 




R 



) 



Ro 



R G annular disc 



for the annular disc in each of the four samples analysed. 

8 Conclusions 

It has been shown, at the level of statistical certainty, that: 

• the optical disc of spiral galaxies consists of two quite distinct dynamical 
regions comprising an interior central sub-disc and a surrounding annular 
disc with an objectively defined dynamical transition boundary separat- 
ing the two regions; 

• that a scaling law, similar to the classical TuUy-Fisher scaling law, applies 
on the interior boundary of the annular disc. 

Furthermore, the analysis has provided a refinement of the Tully-Fisher scal- 
ing law, involving an additional parameter, which renders it insensitive to the 
precise definition of (optical) rotation velocity employed. 



In the following, we describe the analyses corresponding to that of ^ applied 
to the samples of Mathewson & Ford, and Courteaux. For each sample, we use 
two definitions of Vrot- 



• Vrot = Vopt defined as the rotational velocity estimated at Rss- 

For each sample, we find that the use of Vopt in the classical Tully-Fisher scaling 
law gives very significantly worse results that the use of Vauthor- However, we 
then find that performance of the enhanced Tully-Fisher law (with a log Vrot) is 
insensitive to the precise definition used of Vnt with, for example, Vopt being 
comparable to Vauthor for overall effectiveness. 



A The samples of MF and SC and a log Vrot 



• Vr 
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A.l Sample of MFB using Vrot = Vhist 

It is clear from their paper that MFB estimated Vrot for their Tully-Fisher work 
using a version of the histogram method, but it is not clear, either from their 
words or their data, exactly how they implimented the method in practice. In 
the following, we calibrate the basic Tully-Fisher law using MFB's estimates for 
Vrot and absolute magnitudes, Mh say, calculated from Hubble distances with 
H = 85km/ sec /kpc. 

The classical Tully-Fisher scaling law with Vrot = Vhist 

We find, for the classical Tully-Fisher model, 

Mh^Mtf = (-6.72 ± 0.21) + (-6.54 ±0.10) log V,ot; 

(n = 837, tgrad = -(>9) (15) 

Adj r2 = 85%, RSE = 0.45 

after removing 26 of the most extreme outliers. The quoted i-statistic refers to 
the gradient estimate. We note that the gradient lies well within the accepted 
envelope. 

The enhanced Tully-Fisher scaling law with Vrot = Vhist 

We find, for the enhanced model, using exactly the same reduced data set, 

Mh-Mtf = (-7.47 ± 0.23) + [(-6.28 ±0.10) + (0.22 ±0.03) a] log Kot; 

(n = 837, tgrad = -e2, t^ = 6.5) (16) 

Adj = 86%, RSE = 0.44 

so that, according to the ^-statistic for the gradient of the oc log Vrot component 
(tec = 6.5), this component appears to be present in the model in a highly signif- 
icant way. However, each model explains about 85% of the variation in the data 
so that the inclusion of a does not make much difference to model's predictive 
power in this case. 

A.2 Sample of MFB using Vrot = Kpf 

In the following we show that, for the MFB sample, Vrot = Vopt performs poorly 
compared to MFB's own determination, Vrot = V/,;sf , but that when ol is included 
as a predictor, the enhanced model (using Vrot = Vopt) is directly comparable to 
models (|T5|) and ((T6|) above. 
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The classical Tully-Fisher scaling law with Vrot = Vopt 
We find 

Mh ~ Mtf = (-7.70 ± 0.28) + (-5.98 ± 0.13) log Vopf, 
(n = 837, i^,«d = -47) 

Adj = 73%, RSE = 0.60 

after removing twelve of the most extreme outliers. The quoted f-statistic refers 
to the gradient estimate. Note that this model explains only about 73% of the 
variation in the data compared to about 85% for model dTSj) above. 

The enhanced Tully-Fisher scaling law with Vyot = Vopt 

Using exactly the same reduced data set as above, we find 

Mh ~ Mtf = (-8.88 ± 0.22) + [(-5.73 ± 0.10) + (0.78 ± 0.03) a] log Kpf; 

(n= 837, tg,,a = -57, t^ = 26) (17) 
Adj = 85%, RSE = 0.45. 

It is clear that the a. log Vopt component (ta = 26) is very powerfully present in 
the model. Furthermore, this model now explains about 85% of the variation in 
the data (up from 73%) and so is directly comparable with models (|T5l) and (|T6l) 
which both use MFB's own Vrot = Vhist determinations. 

A.3 Sample of Dale et al using Vrot = Vhist 

The primary conclusion from the foregoing analysis of MFB data is that the in- 
clusion of the predictor a In Vrot significantly enhances the performance of the 
Tully-Fisher scaling law. However, it transpires that the degree of this enhance- 
ment depends on the precise definition of Vrot employed, as we now demon- 
strate via the analysis of the DGHU data. 

The classical Tully-Fisher scaling law with Vrot = ^hist 

If we calibrate the basic Tully-Fisher law using Dale's own estimates (0.5 {VgQ% — Viq 
for Vrot and absolute magnitudes calculated from Hubble distances with H = 
85km/ sec /kpc we find, for the TF classical model, 

Mtf = (-7.18 ± 0.28) + (-6.48 ±0.13) log V.ot; (n = 478, t = -52) 
Adj R^ = 85%, RSE = 0.40 
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after removing eighteen of the most extreme outliers. The quoted f-statistic 
refers to the gradient estimate. We note that the gradient lies well within the 
accepted envelope. 

The enhanced TuUy-Fisher scaling law with Vrot = Vhist 

We find, for the enhanced model, using exactly the same reduced data set: 

Mtf = (-7.39 ±0.31) + [(-6.41 ±0.14) + (0.07 ±0.04) a] log 
(n = 478, tgrad = -47, = 1.5) 
AdjR^ = 85%, RSE = 0.40 

so that, according to the i-statistic for the gradient of the a log Vrot component 
(fa = 1.5), this component appears to be not significantly present in the model. 

A.4 Sample of Dale using Vrot = Kpt 

To emphasize the point that it is Dale's definition of Vrot which makes oc log Vrot 
insignificant, rather than some property of the sample, we repeat the analysis 
of DGHU data, but using Vrot = Vopt- 

The classical TuUy-Fisher scaling law with Vrot = Vopt 

We find, for the classical TF model using Vopt — ^(^83)/ 

Mtf = (-7.54 ± 0.38) + (-6.22 ± 0.17) log Vopt) (n = 477, t = -37) 
Adj = 74%, RSE = 0.53 

after removing eighteen of the most extreme outliers. The quoted i-statistic 
refers to the gradient estimate. We note that this model is considerably less 
effective in explaining the data than when Vrot = Vj^ist is used. 

The enhanced TuUy-Fisher scaling law with Vrot — Vopt 
For the enhanced model, we find 

Mtf = (-8.40 ±0.32) + [(-6.01 ±0.14) + (0.61 ±0.04) a] log Vopt; 

(n = 478, f^,,d = -43, t^ = 15) 

Adj R^ ^82%, RSE =0.44 

so that, according to the f-statistic for the gradient of the a log Vrot component 
(fa = 15), this component appears to be present in the model in an extremely 
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significant way. Furthermore, the inclusion of the predictor a log Yyoi has suf- 
ficiently enhanced the predictive power of the model so that it is comparable 
with the model which uses Vrot = Vhist- 

A.5 Sample of MF using Vrot = yhist 

The classical TuUy-Fisher scaling law with Vrot = Vhist 

If we calibrate the basic TuUy-Fisher law using MP's own estimates for Vrot arid 
absolute magnitudes calculated from Hubble distances with H — 85km/ sec /kpc 
we find, for the classical TF model, 

Mtf = (-9.23 ± 0.21) + (-5.50 ± 0.10) log Vrot} (n = 1057, f = -58](18) 
AdiR^^76%, RS£ = 0.46 

after removing 28 of the most extreme outliers. The quoted f-statistic refers to 
the gradient estimate. We note that the gradient is on the low side of expected 
values. 

The enhanced Tully-Fisher scaling law with Vrot = V/j^st 

We find, for the enhanced model, using exactly the same reduced data set, 

Mtf = (-9.88 ±0.23) + [(-5.27 ±0.10) + (0.21 ±0.03) a] log Vrot; 
(n = 1057, tgrad = -53, fa = 6.6) 
AdjR^ ^ 77%, RSE = 0.46 

so that, according to the f -statistic for the gradient of the a. log Vrot component 
{tpc — 6.6), this component appears to be present in the model in a highly signif- 
icant way. 

We see that the inclusion of the predictor a In Vrot gives a significant, but small, 
improvement in performance. 

A.6 Sample of MF using Vrot = Kpt 

The classical Tully-Fisher scaling law with Vrot = Vopt 

If we calibrate the basic Tully-Fisher law using Vopt = ^(^83) we find, for the 
classical TF model, 

Mtf = (-9.83 ± 0.25) + (-5.11 ± 0.11) log Vopt) (n = 1053, f = -47) 
AdjR^ = 67%, RSE = 0.54 
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after removing 32 of the most extreme outliers. We see that the model explains 
about 67% of the variation in the data, compared with about 76% for model ((TSl) . 
It is thus very much less effective. 



The enhanced TuUy-Fisher scaling law with Vrot = Vgpt 

However, if we now consider the enhanced model, we find 

Mtf = (-10.49 ± 0.20) + [(-5.04 ±0.09) + (0.61 ±0.03) a] log Vopf; 

(n = 1053, tgyad = -56, t^, = 23) 

Adj = 78%, RSE = 0.44 

so that, according to the i-statistic for the gradient of the a log Vyot component 
{toi_ = 23), this component is very powerfully present in the model. Furthermore, 
the model now explains about 78% of the variation in the data and is therefore 
significantly better than model (|T8l) which uses MP's own estimate of Vrot- 



A.7 Sample of SC using Vrot = Vmodel = Vmax 

Courteaux's paper is primarily a study of various methods of estimating Vrot for 
TuUy-Fisher studies and his estimator Vnax is derived from as the peak velocity 
achieved by a particular phenomenological model. 



The classical TuUy-Fisher scaling law with Vrot = Vmax 

If we calibrate the basic Tully-Fisher law using Courteaux's Vmax and absolute 
magnitudes calculated from Hubble distances with H = 85km/ sec /kpc we find, 
for the classical TF model, 

Mtf = (-5.53 ± 0.50) ± (-6.66 ± 0.22) log Vrot; {n = 254, t = -30) 
Adj = 78%, RSE = 0.39 

after removing twenty of the most extreme outliers. The quoted i-statistic refers 
to the gradient estimate. We note that the gradient lies well within the accepted 
envelope. 

The enhanced Tully-Fisher scaling law with Vrot = Vmax 

We find, for the enhanced model, using exactly the same data set 

Mtf = (-6.85 ± 0.52) ± [(-6.17 ±0.22) ± (0.36 ±0.06) a] log V,ot; 
(n = 254, tgrad = -28, t^ = 5.9) 

Adj R^ = 81%, RSE = 0.36 
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so that, according to the f-statistic for the gradient of the a log Vrot component 
(ta = 6.3), this component appears to be present in the model in a highly signif- 
icant way. 

A.8 Sample of SC using Vrot = Vmodei = V2.2 

Courteaux's paper is primarily a study of various methods of estimating Vrot for 
TuUy-Fisher studies and his estimator V2.2 is the model velocity at 2.15 x scale 
lengths. 

The classical TuUy-Fisher scaling law with Vrot = V2.2 

If we calibrate the basic TuUy-Fisher law using Courteaux's V2.2 and absolute 
magnitudes calculated from Hubble distances with H — 85km/ sec /kpc we find, 
for the classical TF model, 

Mtf = (-6.43 ± 0.41) + (-6.33 ± 0.18) log Vrot) (« = 262, t = -34) 
Adj = 82%, RSE = 0.36 

after removing eighteen of the most extreme outliers. The quoted i-statistic 
refers to the gradient estimate. We note that the gradient lies well within the 
accepted envelope. 

The enhanced Tully-Fisher scaling law with Vrot = ^2.2 

We find, for the enhanced model, using exactly the same data set 

Mtf = (-6.83 ± 0.47) + [(-6.18 ± 0.20) + (0.10 ± 0.06) a] log Vrot) 
(n = 262, tgrad^-^h tec ^1.8) 
Adj R^ = 82%, RSE =0.36 

so that, according to the ^-statistic for the gradient of the a log Vrot component 
{ta = 6.3), this component appears to be present in the model in a highly signif- 
icant way. 

A.9 Sample of SC using Vrot = Kpt 

The classical Tully-Fisher scaling law with Vrot = Vopt 

If we calibrate the basic Tully-Fisher law using Courteaux's Vopt and absolute 
magnitudes calculated from Hubble distances with H = 85km/ sec /kpc we find. 
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for the classical TF model. 



Mtf = (-6.59 ± 0.64) + (-6.12 ± 0.28) log Vopf, {n = 254, t = -22) 
AdiR^ = 65%, RSE = 0A9 

after removing 21 of the most extreme outliers. The quoted i-statistic refers to 
the gradient estimate. We note that the gradient lies well within the accepted 
envelope. 

The enhanced Tully-Fisher scaling law with Vrot — Vopt 
For the enhanced model, we find 

Mtf = (-7.84 ± 0.51) + [(-5.76 ± 0.22) + (0.75 ± 0.06) oc] log Vopt) 
(n = 254, f« = 13) 

AdjR^^79%, RS£ = 0.38 

so that, according to the i-statistic for the gradient of the a log Vrot component 
{ta — 134), this component is present in the model in an extremely significant 
way. 
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